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Abstract 



Designing an m-bar linkage with a maximal number 
of assembly modes is important in robot kinematics, 
and has further applications in structural biology and 
computational geometry. A related question concerns 
the number of assembly modes of rigid mechanisms as 
a function of their nodes n, which is uniquely defined 
given m. Rigid ll-bar linkages, where n = 7, are 
the simplest planar linkages for which these questions 
were still open. It will be proven that the maximal 
number of assembly modes of such linkages is exactly 
56. 

The rigidity of a linkage is captured by a poly- 
nomial system derived from distance, or Cayley- 
Menger, matrices. The upper bound on the number 
of assembly modes is obtained as the mixed volume 
of a 5 x 5 system. An ll-bar linkage admitting 56 
configurations is constructed using stochastic opti- 
misation methods. This yields a general lower bound 
of f2(2.3") on the number of assembly modes, slightly 
improving the current record of 17(2.289"), while the 
best known upper bound is roughly 4™. 

Our methods are straightforward and have been 
implemented in Maple. They are described in general 
terms illustrating the fact that they can be readily 
extended to other planar or spatial linkages. 
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1 Introduction 

Rigid mechanisms (or linkages) constitute an old but 
still very active area of research in mechanism and 
linkage theory, e.g. [Col02l IFL951 IWH07al IWH07bl 
IWun77j as well as computational geometry and struc- 
tural bioinfor matics, e.g. [BS041 IEM99| IJRKT01| 
IHav98llTD99] . 

A given linkage may be represented by a graph 
with edge set E of lengths kj G R+, for (i, j) e E. 
An assembly mode, or Euclidean embedding, in M d 
is a mapping of its vertices to a set of points in K d , 
such that hj equals the Euclidean distance between 
the images of the i-th and j-th vertices, for € 
E. Euclidean embeddings impose no requirements on 
whether the edges cross each other or not. A linkage 
is (generically) rigid in M. d if and only if, for generic 
edge lengths it can be embedded in M. d in a finite 
number of ways, modulo rigid motions. A graph is 
minimally rigid if and only if it is no longer rigid once 
any edge is removed. 

Let us focus on planar linkages and the associated 
graphs. A graph is called Laman if and only if the 
number of edges is In — 3, where n is the number 
of vertices and, additionally, all of its vertex-induced 
subgraphs on 3 < k < n vertices have < 2k — 3 edges. 
This is essentially the Griibler-Kutzbach-Chebychev 
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Figure 1: An 11-bar linkage. 

formula on the degrees of freedom for mechanical 
linkages, e.g. |Ang89| . It is a fundamental theorem 
that the class of Laman graphs coincides with the 
generically minimally rigid planar linkages |Lam70] . 
The problem studied in this paper is to compute the 
number of distinct planar assembly modes of rigid 
mechanisms, up to rigid motions. In particular, the 
maximal number of assembly modes of an 11-bar link- 
age (Figure [1]) will be presented for the first time. 

1.1 Existing work. 

The algebraic approach is to define a well-constrained 
polynomial system expressing the length constraints, 
such that the system's real solutions correspond pre- 
cisely to the different assembly modes. When defin- 
ing a straightforward system such as ((T), where the 
unknowns are the nodes' coordinates, all nontrivial 
equations are quadratic. For planar linkages, there 
are In — 4 equations hence, by applying the classical 
Bezout bound on the number of common roots, we 
obtain 4™~ 2 . It is indicative of the hardness of the 
problem that efforts to substantially improve these 
bounds have failed. 

Today, the best general upper bound is roughly 
4"- 2 /V 7r O - 2). This was obtained using determi- 
nantal varieties defined by distance matrices BS04 . 
Straightforward application of mixed volumes (dis- 
cussed in Section [3) yields an upper bound of 4 n ~ 2 

[sTToj . 

The best general lower bounds are 24 L(" — 2 )/ 4 J ~ 



Figure 2: The Desargues linkage or, equivalently, a 
planar parallel robot. 

2.21" and 2 ■ i2L(«-3)/3j ~ 2.289"/6, obtained 
by a caterpillar and a fan construction, respec- 
tively [BS04J . Both bounds are based on the De- 
sargues graph (Figure |2J| , which admits 24 assembly 
modes. 

In applications, it is crucial to know the number of 
assembly modes for specific (small) values of n. The 
most important result in this direction was to show 
that the Desargues linkage, also known as the planar 
parallel robot, admits precisely 24 assembly modes 
in the plane [BS041 IHun83l IGSR921 ILM94] , More- 
over, the if 3,3, or 9-bar, linkage admits exactly 16 
assembly modes in the plane WHO 7a . This paper 
demonstrated such a linkage through an elaborate 
construction which was not required in our case. 

All existing bounds for n < 10 are found in Table [1] 
from |ETV09j . 

1.2 Our contribution. 

To upper bound the number of planar rigid mecha- 
nisms, we explore adequate polynomial systems lead- 
ing to tight root bounds. Modeling physical systems 
by appropriate polynomial systems is a deep and hard 
question, with a wide range of applications in differ- 
ent fields. 

We employ two powerful algebraic tools for defin- 
ing polynomial systems and for counting the sys- 

x This slightly corrects the exponent of the original state- 
ment. 



2 



11 = 


3 


4 


5 


6 


7 


8 


9 


10 


upper 


2 


4 


8 


24 


64 


128 


512 


2048 


lower 


2 


4 


8 


24 


48 


96 


288 


576 



Table 1: Bounds for the number of embeddings of rigid graphs with n < 10 [ETV09] . 



tern's common roots. For the former, we use dis- 
tance matrices, also known as Cayley-Menger matri- 
ces, which contain all known and unknown distances 
between the graph's nodes. The signs of the ma- 
trix minors capture rigidity and embeddability in Eu- 
clidean spaces, as described in Section |3l Our results 
indicate that such matrices are advantageous to us- 
ing the coordinates' formulation, such as system (JXJ) , 
when constructing the polynomial system. 

Our second tool is the mixed volume of a well- 
constrained polynomial system, which exploits the 
sparseness of the equations. The mixed volume 
bounds the number of common roots, by Bernstein's 
Theorem [TJ as described in Section [3] This bound 
is never larger than Bezout's, and is typically much 
tighter for systems whose equations do not contain all 
possible terms for a given total degree. It turns out 
that such are the systems encountered here and could 
be of wider interest in similar enumeration problems. 

Moreover, mixed volume is a bound so it only needs 
to consider which are the nonzero coefficients, with- 
out considering specific values. This is in sharp con- 
trast to solving methods often used for root counting, 
where one assigns random values to the coefficients. 

We thus obtain a tight bound for linkages with 
n = 7 nodes by appropriately formulating the poly- 
nomial system based on distance matrices. We com- 
pute the mixed volume of a 5 x 5 polynomial system, 
which equals 56. This is tight since we demonstrate a 
construction with as many assembly modes, see Fig- 
ure [5] We also show how our stochastic optimisation 
methods are generalizable for studying further pla- 
nar or spatial linkages. Our construction also yields 
a general lower bound of Q(2.3 n ) on the number of 
configurations, thus improving the available bound of 
17(2.289"), whereas the best available upper bound is 
0(4"/Vi). 

The rest of the paper is structured as follows. 



Section [2] presents rigid mechanisms and how our 
study focuses on a specific 11-bar linkage. Section [3] 
presents our algebraic tools and obtains the upper 
bound of 56 for 7 nodes and 11 bars. Section [4] ac- 
tually constructs a linkage with 56 assembly modes. 
We conclude with open questions. 

2 Rigid linkages 

This section studies the number of assembly modes 
of planar rigid linkages or, equivalently, rigid graphs 
and their number of embeddings in R 2 . 

Such linkages admit inductive constructions that 
begin with a triangle, followed by a sequence of so- 
called Henneberg steps. There are two types of such 
steps, each adding one new vertex and increasing the 
total number of edges by 2. A graph is Laman, or the 
associated mechanism is rigid, if and only if it can 
be constructed by a sequence of the corresponding 
Henneberg steps. By exploiting this fact, all rigid 
graphs in M. 2 , but also in K 3 , were constructed using 
the Henneberg steps in [ETV09] , and classified up to 
graph isomorhism. 

Let us consider the Henneberg steps defining 
Laman graphs, each adding a new vertex and a total 
number of two edges. A Henneberg- 1 (or Hi) step 
connects the new vertex to two existing vertices. A 
Hennenerg-2 (or H2) step connects the new vertex 
to 3 existing vertices having at least one edge among 
them, and this edge is removed. Both steps are illus- 
trated in Figure [3] We represent each Laman graph 
by AS4 . . . , s„, where € {1, 2}; this is known as its 
Henneberg sequence. Note that this sequence is by 
no means unique. A Laman graph is called H\ if and 
only if it can be constructed using only H\ steps; it 
is called H2 otherwise. 

Since two circles intersect generically in two points, 
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Figure 3: The planar Henneberg steps; the bottom 
vertex is new. 

a Hi step at most doubles the number of assembly 
modes and this is tight, generically. In Table Q] most 
lower bounds follow from this fact. The lower bound 
for n — 9 follows from the Desargues fan [BS04 . 

For a Laman graph on 6 vertices, a tight uppper 
bound of 24 follows by examining the only three pos- 
sibilities: the graph is either Hi, it is K 3 , 3 , or is the 
Desargues graph. Now, Hi graphs on 6 vertices have 
exactly 16 assembly modes. The K 3 ^ 3 graph has pre- 
cisely 16 assembly modes [WH07aj . a fact first con- 
jectured in |Wun77j . The Desargues graph has pre- 
cisely 24 assembly modes: the upper bound was first 
shown in |Hun83] and proven more explicitly, along 
with the lower bound, in |GSR92j . This fact was in- 
dependently rediscovered in |BS04) and proven via 
coupler curves. 

In [KTV09) further upper bounds were obtained 
from the following family of simple systems, employed 
in order to express graph embeddability in R 2 . Let 
Xi , yi denote the coordinates of the i-th vertex, where 
the Uj are the given lengths: 

( xi = Oi, yi = h, i = l,2, ai,bi€R 

I fa - x 3 ) 2 + ( Vi - yj) 2 = % (i, j) eE- {(1, 2)} 

(1) 

Vertex (xi,yi) was fixed to (0, 0), to discard transla- 
tions, and (#2,2/2) t° (1)0) to remove rotations and 
scaling, assuming without loss of generality that the 
corresponding edge exists. Then the mixed volume of 
the system gave the upper bounds in the table which, 
unfortunately, are loose, even after removing roots at 
infinity [ETV09 . Below we overcome this limitation 
by considering polynomial systems derived from dis- 
tance matrices. 

We now study the case n = 7. If a Hi step is ap- 
plied to any linkage with n = 6, the resulting linkage 
with n = 7 admits exactly 48 assembly modes, which 
is reflected to the corresponding lower bound in Ta- 



bic [U To maximize the number of assembly modes, 
we shall apply a Hi step to any graph with n = 6. 
By checking graph isomorphisms and taking into ac- 
count the various symmetries, it was shown |ETV09| 
that there are only 3 relevant graphs to be consid- 
ered. These are obtained by a H 2 step applied to the 
Desargues graph as follows, where notation refers to 
Figured 

We remove edge (4,5), and add edges (4, 7), (5, 7), 
and 

• (1,7), or 

• (3,7), or 

• (6,7). 

The first case corresponds to the topology that 
shall be studied extensively in the sequel, since it 
leads to the maximum number of 56 assembly modes. 
The linkage is shown in Figure refllbarLinkage and 
in the next section we will proove that 56 is the max- 
imum number of assembly modes possible for such 
a linkage. The other two linkages admit at most 44 
and 48 assembly modes, respectively, hence they are 
not studied any further. These upper bounds are ob- 
tained as mixed volumes of polynomial systems, as 
explained in Section [3] 

3 Algebraic methods for the 
upper bound 

This section discusses polynomial systems, describes 
distance matrices, and derives an upper bound on the 
number of assembly modes. 

3.1 Mixed volume 

We first discuss multivariate polynomial systems 
and introduce sparse elimination theory in order to 
exploit their structure and sparseness; for details, 
see [CLO05, Emi94 . Classical elimination theory 
characterizes every polynomial by its total degree. 
For a well-constrained system of polynomial equa- 
tions, the classical Bczout bound on the number of 
isolated roots equals the product of the polynomials' 
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total degrees. One disadvantage of this bound is that 
it counts projective roots and hence increases when 
there are roots at projective infinity. 

In sparse (or toric) elimination theory, a polyno- 
mial is characterized by its support. Given a poly- 
nomial / in n variables, its support is the set of ex- 
ponents in N n corresponding to nonzero terms (or 
monomials) . The Newton polytope of / is the convex 
hull of its support and lies in M. n . Consider (Newton) 
polytopes p C R™ and parameters Ai G l,Aj > 0, 
for i = 1, . . . , n. We denote by XiPi the correspond- 
ing scalar multiple of Pi. Consider the Minkowski 

sum of the scaled polytopes AiPi H h A n P„ € R n ; 

its Euclidean volume is a homogeneous polynomial of 
degree n in the A^. The coefficient of Ai - • • A„ is the 
mixed volume of Pi, . . . , P n . If Pi = ■ ■ ■ = P n , then 
the mixed volume is n! times the Euclidean volume 

of Pi. 

We now focus on the topological torus C* = C— {0} 
in order to state Bernstein's root bound in terms of 
mixed volume. 

Theorem 1 [Bcr75 Let fi = ■ ■ ■ = f n = be a poly- 
nomial system in n variables with real coefficients, 
where the fi have fixed supports. 

The number of isolated common solutions in (C*) n 
is bounded above by the mixed volume of the Newton 
polytopes of the fi . This bound is tight for a generic 
choice of coefficients of the fi 's. 

One alternative to using general bounds such as 
mixed volume is to manipulate a system, with coef- 
ficients chosen randomly, so as to bound the number 
of common real roots, e.g. by means of Grobner bases 
or homotopy continuation. We actually used Grobner 
bases for all systems studied in this paper, where we 
have used random coefficients. In particular, for the 
linkage analyzed below, we computed the total-degree 
Grobner basis of system (jTJ) with random distances 
in Maple, and obtained the Hilbert polynomial of the 
corresponding ideal. 

This is a univariate polynomial where, by setting 
its variable to 1, one deduces an upper bound on the 
number of complex common roots. This was indeed 
56, which coincides with mixed volume. However, 
this only offers an indication on the number of roots, 



not an actual bound, since it depends on the choice 
of coefficients. 

The advantage of mixed volume is that it treats 
entire classes of systems defined by their nonzero 
terms, without considering specific coefficient val- 
ues. Lastly, a tight mixed volume implies that one 
can solve the system efficiently either by sparse re- 
sultants, e.g. [CEOOj . or by sparse homotopies, e.g. 
[Ver99] . 

3.2 Distance geometry 

The theory of distance geometry has been well devel- 
oped, e.g. |Blu701 Sch35] , with several applications 
especially in structural bioinformatics, e.g. EM99. 
IHav98j . 

Definition 1 Given an n-vertex graph, the corre- 
sponding distance or Cayley-Menger matrix is a sym- 
metric (n + 1) x (n + 1) matrix B, indexed by i = 
0, . . . , n, such that B(0, i) = 1 for i > 1, B{i, i) = 
for i > 0, and B{i,j) equals the squared distance be- 
tween vertices i,j, for j > i > 1. 

To be more precise, to obtain the Cayley-Menger 
matrix of the given n vertices, we must multiply B 
by —1/2. But for testing embeddability, the two for- 
mulations are equivalent, as seen in the following 
theorem, due to the work of Cayley and Menger, 
see |Blu70[ ISch35) . Let D(ii,...,i k ) denote the 
(k + 1) X (k + 1) diagonal minor of B indexed by 
rows and columns 0, i\, . . . , where ii < • • • < i k € 

Theorem 2 Suppose we are given a matrix B of the 
form specified in Definition [3 Then B corresponds 
to a (complete) graph embeddable in M. d , if and only 
if 

1. for k = 2, . . . , d + 1 and any {i\, . . . , i k } C 
{1, . . . , n}, it holds {—l) k D{ii, . . . ,ik) > 0, and 

2. rank(B) =d + 2. 

The second condition, due to Cayley, yields a 
(large) number of equalities, which are typically not 
independent. 
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The first condition, due to Monger, yields inequal- 
ities: For k = 2, it expresses the fact that all entries 
must be non-negative. For k = 3, it captures the tri- 
angular inequality. If we apply it, without loss of gen- 
erality, to indices 1,2,3, using Cij = ifj for the given 
entries, then the condition states that D(l, 2, 3) < 0, 
where this quantity equals 

— ('12+^13+^23) ('12+^13—^23) (I12+I23— 113) (^13+^23~^12 

This vanishes precisely when the corresponding 
points are collinear, whereas D(l, 2, 3) < when the 
points define a triangle. Equivalently, D(i,j,k) < 
can be written 1^ + Ijk > hj for all triplets k € 
{l,...,n}. For k = 4 the condition captures the 
tetrangular inequality. 

3.3 Upper bound 

We shall employ the rank condition in order to derive 
equality constraints on the unspecified distances of 
our linkage. Here is the matrix, where the c 2 j = If^ 
correspond to the fixed distances, and Xij are the 
unspecified distances. 









i' 2 


"3 


V 4 


v 5 


!'<> 


V7 




' 


1 


1 


1 


1 


1 


1 


1 


Vl 


1 





Cl2 


Cl3 


C14 


Xl5 


xie 


Cl7 


V2 


1 


C12 





C23 


2^24 


C25 


X26 


2^27 


"3 


1 


Cl3 


C23 





2^34 


2^35 


C36 


237 


U 4 


1 


C14 


2^24 


£34 





X45 


f'46 


C47 


"5 


1 


X15 


C25 


2:35 


2:45 





C56 


C57 


"6 


1 


xw 


X26 


C36 


C46 


C56 





2^67 


v 7 


1 


C17 


X27 


2^37 


C47 


C57 


X67 






By Theorem [21 any 5x5 minor of this matrix must 
vanish, which yields polynomial equations on the 
variables x^. There are (3) = 35 such minors, each 
in 2 to 4 variables. Among these polynomials, no 
4x4 subsystem exists that corresponds to a rigid 
mechanism, as can be verified by checking Laman's 
condition. This is indispensable for the system to 
have a finite number of solutions. 

However, it is possible to find certain 5x5 sub- 
systems whose subgraph is Laman and, moreover, 
uniquely define the configuration of the overall link- 
age (Figure 0]) . One of these systems has 4 bivariate 



equations and one trivariate equation, and is defined 
by taking the following diagonal minors: 

' D(4, 5,6, 7)(C46,C47,C 56 ,C 57 ,X45,2; 67 ) = 
D(l, 4, 6, 7)(ci4, C17, C 46 , C47, 27i 6 , x 67 ) = 

< D(l,4, 5,7)(ci4,ci 7 ,C47,c 57 ,a;i5,X45) = (2) 
D(l, 2, 3, 5)(ci2, ci 3 , c 25 , c 23 , xi 5 , x 35 ) = 
)• [ ^(1,3, 5, 6)(ci 3 ,C36,c 56 ,a;i5, ^16,2:35) = 

defines a system of 3 quadratic and two cubic equa- 
tions in X15, ccie, 2:35, X45, £67; the cubics are the first 
and last polynomials. Now, this system's mixed vol- 
ume turns out to be 56, computed using the software 
from |EC95j 0. This bounds the number of (com- 
plex) common system's solutions, hence the number 
of (real) assembly modes. 

Notice that this bound does not take into ac- 
count solutions with zero coordinates, in other words 
some zero length. However, a linkage has assembly 
modes with some zero length only when the input 
bar lengths form a singular set, in the sense that 
they would satisfy a non-generic algebraic depen- 
dency. For example, by letting some input distance 
be exactly 0, some 11-bar linkage may theoretically 
have infinitely many configurations. However, gener- 
ically, it is impossible to have such a linkage. 

Thus, the mixed volume bound guarantees that a 
true 11-bar linkage can have at most 56 assembly 
modes. 

4 Lower bound 

To prove that 56 is a lower bound on the number 
of configurations, it is sufficient to exhibit a linkage 
with 56 configurations. Unfortunately, we have 11 
design parameters. Thus our search space is homco- 
morphic to R^ 1 . Even if we reduce the search space 
to the integers between 1 and 100 for each parameter, 
an exhaustive search would lead us to consider 10 22 
linkages. 

Another approach was used in [LM94J to find a pla- 
nar parallel robot with 12 real solutions. The authors 
found a singular positions with many real solutions, 

2 http://www. di.uoa.gr/ emiris/ 
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4.1.1 Evaluation of the objective function 




Figure 4: When the edge lengths [15], [16], [35], [45] 
and [67] are specified as a solution of System ^ , the 
11-bar linkage has only one assembly mode. 



and use deformations to get the desired result. The 
advantage of this approach is that the space of singu- 
lar configurations has a lower dimension than the full 
design space. However, in our case, the dimension 
was still too high and this approach could not allow 
us to conclude. 

4.1 Stochastic Methods 

Random sampling is a common approach to search 
in high dimensional space is random sampling. The 
system at hand being homogeneous, we can specify 
the last coordinate Zio to 100 without restriction of 
generality. Moreover, in the following, we restrict our 
domain to integer coordinates. 

Our goal is to maximize the function: 



N : 



(lo, Z9) 



f 0, if System @ has infinitely 
1 many solutions 



The cost of evaluating the objective function is the 
bottleneck of the stochastic methods. The function N 
needs to be evaluated several hundreds of time before 
converging toward a point where it is maximal. 

In our case, evaluating N means solving a sys- 
tem of polynomial equations with finitely many so- 
lutions, and counting the number of its real posi- 
tive roots. Experimentally, we use the Maple func- 
tion RootFinding [Isolate] . This function com- 
putes isolating boxes around the real solutions of the 
system to solve. The underlying algorithm uses Grb- 
ner bases and Rational Univariate Representation. 

The choice of the system modeling the 11-bar ma- 
nipulator is critical. A naive modelling yields a sys- 
tem with 14 variables, one for each of the coordinates 
of the 7 vertices, and 11 equations, one for each link 
in the manipulator. In Maple, solving this system 
is roughly 10 times slower than solving System ([2|). 
That can be explained partly by the fact that, gener- 
ically, the time complexity to solve a 0-dimensional 
system is roughly exponential in the number of vari- 
ables, e.g. jCLO05llLai0l] . 



4.1.2 Simulated annealing 

The Monte-Carlo simulation and its simulated an- 
nealing variant [KGJV83] are the most well spread 
stochastic optimisation methods. These methods 
have already been used in robotics for path planning 
in [KSL096 . The convergence of these methods has 
been well studied and the simulated annealing sim- 
ulations has been successfully used in different field 
such as biology and chemistry. We implemented this 
method to maximize the function N. 



I Number of solutions of System The Monte-Carlo simulation depends on a parame- 



otherwise 



A first idea is to sample random points using a 
Gaussian law centered on an arbitrary initial point. 
This naive approach did not allow us to find a linkage 
with 56 configurations, but we could finally find such 
a point by using more sophisticated methods. 



ter T called temperature. In the simulated annealing 
variant, the temperature T is decreased according to 
a specific schedule at each step. The optimal way 
to decrease the temperature depends on each prob- 
lem. In our case, we chose arbitrarily a linear cooling 
schedule ( [SH87|lNA"98] ). The corresponding pseudo- 
code is summarized in Algorithm [T] 
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la = 180 h = 70 l 2 = 200 l 3 = 205 
h = 210 l 5 = 205 l 6 = 80 Z 7 = 200 
l 8 = 70 Z 9 = 200 l 10 = 100 

has exactly 56 assembly modes. Some of its assem- 
bly modes are shown in Figure [5j 

Our simulations show that the set of 11-bar link- 
ages with 56 assembly modes has a non-zero volume. 
With the direct random sampling, we did not find 
any linkage with more than 48 assembly modes. This 
indicates that linkages with bar lengths between 
and 300 have a very small probability to have 56 as- 
sembly modes. Our experiments show that the Cross 
Entropy method is a general approach well suited to 
generate such configurations. Moreover, it can be 
used easily for other mechanisms, such as larger pla- 
nar or spatial linkages for example. A variant of this 
method can also be used to compute an estimation 
of the probability to find a linkage with a maximal 
number of assembly mode. 



4.1.3 Cross entropy method 

A new method was introduced recently by Rubin- 
stein in |Rub97j for the simulation of rare events. 
This method is especially well-suited for combinato- 
rial and continuous optimization ([Rub99 ). The idea 
of this approach consists in minimizing the distance 
between specific probability laws appearing during 
the computation (see [DBKMR05] for more details). 
For our problem, we use the scheme developed for 
continuous multi-extremal function optimization pre- 
sented in |KPR06j . The funct ion N is not continu- 
ous, but behaves smoothly: generically, if we modify 
the lengths of a given linkage by small enough values, 
the number of configurations is increased or decreased 
only by 2. For our problem it yields Algorithm [2j 



4.2 Results 

In Table [2 we compare the results of three stochas- 
tic methods used to optimize the function N with 
respect to the lengths of the 11-bar linkage. 

The first column shows the result of the direct ran- 
dom sampling. Each line correspond to 600 evalua- 
tion of N on points chosen according to a Gaussian 
law of center 100 and standard deviation 100. This 
operation has been run 10 times, yielding a total of 
6000 evaluations of N on random points. However, 
this did not allow us to find a linkage with 56 config- 
urations. 

The second approach is a Simulated Annealing sim- 
ulation. We ran it 10 times, limiting the number of 
evaluations of N to 600. One of our simulation re- 
turned a set of link lengths yielding a linkage with 
56. 

Finally, the last column reports the results of the 
Cross Entropy Method. Each line correspond to a 
simulation stopped after 600 evaluations of N. This 
approach was the most efficient. Four simulations out 
of ten returned a linkage with 56 configurations. 

The results of our stochastic methods are summa- 
rized in Table [2] In particular, we found that the 
manipulator with the following design parameters 



Direct Sampling Simulated Anealing Cross Entropy 



44 (572) 


52 (17) 


52 (199) 


42 (196) 


54 (247) 


54 (132) 


48 (27) 


48 (362) 


52 (186) 


44 (200) 


52 (14) 


54 (130) 


42 (200) 


54 (547) 


56 (497) 


44 (424) 


54 (315) 


56 (328) 


46 (48) 


56 (425) 


56 (454) 


42 (170) 


50 (585) 


54 (375) 


42 (18) 


54 (26) 


56 (552) 


46 (366) 


52 (474) 


56 (355) 



Table 2: Results of different stochastic optimisation 
algorithms. Each line corresponds to a simulation 
stopped after 600 evaluations of N. The results is of 
the form n(m) where n is the maximal reached value 
of N and m is the number of evaluations done before 
reaching it. 
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Algorithm 1 Simulated annealing with linear cool- 
ing to maximize the number of solutions of System 



(2) 



T ^4 

max Step <— 1000 

U i — [10, ...,10] 

point <- {l = 100, h = 100} 
value 4— N(lo, lg) 

for n from 1 to maxStep while value < 56 do 

newPoint 4— GaussianN eighbour (point, a) 
newValue 4— N(Iq, lg) 
treshold 4— Uniform random value in [0,1] 
newValue > value 



if < 



or 



then 



treshold < e 
point <— newPoint 
value <— newValue 
end if 

T <— To(l — nj max Step) 
end for 
return point 



Algorithm 2 Cross Entropy Method to maximize 
the number of solutions of System @ 

maxStep <— 600 

Nehte <~ 5 

N <- 20 
a <- 0.5 

[100, ...,100] 
<r^- [100, ...,100] 
7^0 
max <— 

for n from 1 to maxStep while max < 56 do 
for i from 1 to N do 

Xi <— GaussianN eighbour((l, sigma) 
end for 

prevmax 4— max 
for i from 1 to N do 
V, <- AT(Xi) 
if Vi > max then 
max 4— V, 
point <— Xj 
end if 
end for; 

Sort X and V from the largest V$ to the small- 
est 

7 <- Vjv aMt . 

r+new , 1 \~^N e i ite ^ 



(T 



new 



^£f""-(x,-£) 2 



A, 



/i <- a/I ne "' + (1 - a) /2 
(7 4- a<7 neuj + (1 - a)<7 
end for 
return point 
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5 Improvement on the general 
lower bound 

This specific linkage allows us to improve slightly the 
lower bound on the maximal number of assembly 
modes of a planar linkage with n pivot joints. We 
use the same fan construction as in [BS04 . 

We construct a linkage with n nodes composed of 
k sub-mechanisms as follows. Each of the k sub- 
mechanisms is a 11-bar linkage with the topology we 
have considered so it has 56 assembly modes, and 
28 when one triangle is fixed. All sub-mechanisms 
share this triangle so there are no degrees of free- 
dom other than those within each sub-mechanism, 
and the overall mechanism is rigid. The total number 
of assembly modes equals the product of remaining 
assembly modes per sub-mechanism, since these sub- 
mechanisms can be configured independently. Hence 
the overall number is 28 k when one triangle is speci- 
fied, and 2 • 28 k in total. For n pivot joints, we let k 
be the greatest integer < ^j^. This yields the lower 
bound: 

2 • 28 l ^ i 

The constant under the exponent is v^28 ~ 2.3003 
while in the previous lower bound the constant was 
\/Y2 ~ 2.2894. 



6 Further work 

Undoubtedly, the most important and oldest problem 
in rigidity theory is the full combinatorial character- 
ization of rigid graphs in R 3 . We believe that one 
can extend our methods to spatial linkages. We ex- 
pect our approach leads to an algorithmic process for 
obtaining good algebraic representations of the enu- 
meration problem at hand, namely low mixed vol- 
umes for the systems derived from distance matrices. 
One issue is that the number of equations produced is 
quite large, with algebraic dependancies among them. 
The question becomes then to choose the best well- 
constrained system. 

The structures studied in this paper are point- 
and-bar structures; they generalize to body-and-bar, 
where edges can be connected to different points of a 



rigid body. It is known that a body-and-bar structure 
in R d is rigid if and only if the associated graph is the 
edge-disjoint union of C^ 1 ) spanning trees [WT85| . 
Body-and-bar structures are our next object of study. 
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SAGA: Shapes, Algebra, and Geometry. He also ac- 
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